Bogdat etal. BMC Systems Biology 201 3, 7:67 
http://www.biomedcentral.eom/1752-0509/7/67 



Systems Biology 



RESEARCH ARTICLE 



Open Access 



Levels of pro-apoptotic regulator Bad and 
anti-apoptotic regulator Bcl-xL determine the 
type of the apoptotic logic gate 

Marta N Bogdat 1 ' 2 , Beata Hat 1 , Marek Kochahczyk 1 and Tomasz Lipniacki 1 ' 3 * 



Abstract 

Background: Apoptosis is a tightly regulated process: cellular survive-or-die decisions cannot be accidental and 
must be unambiguous. Since the suicide program may be initiated in response to numerous stress stimuli, signals 
transmitted through a number of checkpoints have to be eventually integrated. 

Results: In order to analyze possible mechanisms of the integration of multiple pro-apoptotic signals, we 
constructed a simple model of the Bcl-2 family regulatory module. The module collects upstream signals and 
processes them into life-or-death decisions by employing interactions between proteins from three subgroups of the 
Bcl-2 family: pro-apoptotic multidomain effectors, pro-survival multidomain restrainers, and pro-apoptotic single 
domain BH3-only proteins. Although the model is based on ordinary differential equations (ODEs), it demonstrates 
that the Bcl-2 family module behaves akin to a Boolean logic gate of the type dependent on levels of BH3-only 
proteins (represented by Bad) and restrainers (represented by Bcl-xL). A low level of pro-apoptotic Bad or a high level 
of pro-survival Bcl-xL implies gate AND, which allows for the initiation of apoptosis only when two stress stimuli are 
simultaneously present: the rise of the p53 killer level and dephosphorylation of kinase Akt. In turn, a high level of Bad 
or a low level of Bcl-xL implies gate OR, for which any of these stimuli suffices for apoptosis. 

Conclusions: Our study sheds light on possible signal integration mechanisms in cells, and spans a bridge between 
modeling approaches based on ODEs and on Boolean logic. In the proposed scheme, logic gates switching results 
from the change of relative abundances of interacting proteins in response to signals and involves system bistability. 
Consequently, the regulatory system may process two analogous inputs into a digital survive-or-die decision. 

Keywords: Apoptosis, Cell survival, Signaling pathway, Bcl-2 family, Bistability, Boolean logic, Ordinary differential 
equations 



Background 

Biological background 

Apoptosis is a programmed cell death required for 
removal of infected, damaged or unwanted cells [1]. It 
assists in the development and aging as a homeostatic 
mechanism controlling cell populations in tissues, and 
it constitutes a key immune defense mechanism against 
infected or cancerous cells [1,2]. Disrupted regulation of 
apoptosis contributes to carcinogenesis, spread of infec- 
tion, autoimmunological and neurodegenerative disorders 
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[3-5]. Since there could be many reasons for which a 
cell should be eliminated, there exist numerous pathways 
(major of them discussed below) through which apoptosis 
can be initiated. This opens an interesting question how 
the pro-apoptotic, but also pro-survival, signals converge 
and are integrated before the survive-or-die decision is 
reached. 

In mammalian cells, apoptosis can be induced via 
two classes of apoptotic pathways: extrinsic and intrin- 
sic [6,7]. In both classes, signaling leads to the activation 
of a family of cysteine proteases named caspases which 
serve as executors of the apoptotic process [8,9]. Cas- 
pases are present in virtually every cell in the form of 
inactive precursors called pro-caspases [9]. Each apop- 
totic pathway activates some initiator caspases which, in 
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turn, activate the main executioner, caspase-3 [8]. The 
release of cytochrome c from mitochondria results in 
the formation of apoptosome (containing also Apaf-1 
and caspase-9) and activation of caspase-9 [10], which 
activates caspase-3 triggering the caspase cascade. This 
cascade involves caspase-2, -6, -8 and -10 [11,12] respon- 
sible for the proteolytic dismantling of the apoptotic cell 
[13]. Caspase-3 is responsible for the further release of 
cytochrome c, which ensures that the apoptotic decision is 
irreversible [14,15]. 

The extrinsic apoptotic pathways, also known as 
'death receptor pathways', are initiated through the acti- 
vation of membrane death receptors, including Fas, 
TNFR, DR3 and DR4/DR5, by their respective ligands 
(FasL, TNFa, TWEAL and Trail) [16,17]. Engagement 
of death receptors by their cognate ligands triggers 
the recruitment of different adaptor proteins. Depend- 
ing on the recruited adaptors, either pro-apoptotic sig- 
nals (mediated by caspases-8 and -10) or pro-survival 
(MAPK- and/or NF-KB-mediated) signals are induced. 
Caspase-8 and -10 recruited to death receptors autoac- 
tivate themselves and activate further effector caspases 
(caspase-3 and -7) either by direct processing or by engag- 
ing the intrinsic death pathway [11]. 

The intrinsic apoptotic pathway is also called the mito- 
chondrial pathway because it is associated with the 
disruption of mitochondrial outer membranes and conse- 
quent release of cytochrome c. This process is regulated 
by the Bcl-2 family proteins, which share up to four BH 
(Bcl-2 homology) domains. The proteins can be classified 
as: pro-apoptotic multidomain effectors (Bax and Bak), 
pro-survival multidomain restrainers (including i.a. Bcl-2 
proper, Bc1-xl, Bcl-w, Mcl-1, Al) and pro-apoptotic sin- 
gle domain BH3-only upstream sentinels (comprising Bid, 
Bim, Bad, Puma, Noxa and others) [18,19]. 

Bax and Bak are effectors directly responsible for 
the mitochondrial outer membrane permeabilization via 
either channel formation [20] or opening of voltage- 
dependent anion channels [21]. While these effectors are 
present even in surviving cells, activated p53 can induce 
transcription of their genes and further elevate levels of 
both Bax [22] and Bak [23]. Deletion of either bax or 
bak affects apoptosis only slightly, but deletion of both 
these genes dramatically impairs apoptosis in many tissues 
[24,25]. Despite intense studies, it is still controversial, 
how the level, conformation and activity of these pro- 
apoptotic effectors is regulated. There is a bulk of evidence 
that Bax and Bak (plausibly after initial conformational 
priming by BH3-only proteins [26,27]) can be inhibited by 
pro-survival restrainers [28-30] (either by direct seques- 
tration in the mitochondrial membrane [31] or by active 
retrotranslocation to the cytosol [32,33]). Other studies 
suggest that Bax and Bak could be activated directly by 
some BH3-only proteins (Bid, Bim, Puma) [26,27,34,35], 



however bim~'~ bid~'~ mice develop rather normally [36] . 
Additionally, the release of Bax from pro-survival restrain- 
ers appears to be sufficient to trigger apoptosis [37], 
and these interactions have been characterized in detail 
recently [38,39]. The function of another subset of BH3- 
only proteins (Bad, Bid, Noxa) is proposed to rely mainly 
on displacing (conformationally primed) pro-apoptotic 
effectors from pro-survival restrainers [40,41]. This mech- 
anism is enabled by much higher affinity of BH3-only pro- 
teins to restrainers in comparison to effectors: in cells with 
abundant Bcl-XL:Bad heterodimers, no Bax:Bcl-XL het- 
erodimers are present [31,42]. In normal cells Bad is found 
mostly in the phosphorylated form (at Serll2 and Serl36 
by phosphorylated Akt), in which it preferentially binds to 
the scaffolding protein 14-3-3 (isoform 9) instead of Bcl- 
2 family pro-survival restrainers [43-45]. Although direct 
interaction of Bax and 14-3-3 has been demonstrated and 
overexpression of 14-3-3 inhibits apoptosis, prior cas- 
pase activity is required for cleaving 14-3-3 to release 
Bax [46]. 

Bcl-2 family proteins are activated or inhibited in 
response to numerous stress factors including heat shock, 
y and UV irradiation, nutrient deprivation, viral infection, 
hypoxia and increased intracellular calcium concentration 
[47,48]. In this study, we confine to two sources of exter- 
nal stimuli: 1) DNA damage prompting the activation of 
p53 (which mediates pro-apoptotic signals) [49,50] and 2) 
withdrawal of growth factors (GF) leading to the deactiva- 
tion of Akt (which, when active, mediates anti-apoptotic 
signals) [51-53]. 

In not onco-transformed cells, p53 protein remains 
inactive [54]. In response to DNA damage, p53 is acti- 
vated by phosphorylation at Serl5 and Ser20 (by ATM 
[55,56]), which protects it from rapid degradation [56,57]. 
When phosphorylated, p53 is capable of inducing synthe- 
sis of its own inhibitors: ubiquitin-protein ligase Mdm2 
and serine/threonine phosphatase Wipl [58,59], as well 
as proteins responsible for cell cycle arrest and DNA 
repair [60]. Additional p53 phosphorylation at Ser46 
(by kinase HIPK2 [61,62]) enables p53 to activate expres- 
sion of proteins which mediate apoptosis [63], in par- 
ticular pro-apoptotic Bax and Bak [22,23]. DNA repair 
and apoptotic functions make p53 a primary tumor sup- 
pressor; respectively the p53 gene is the most frequently 
mutated gene in cancers [64,65]. 

In healthy cells but also frequently in cancer cells, 
Akt, in contrast to p53, maintains its (at least partial) 
activity and suppresses apoptotic signals by phosphory- 
lating and thereby inhibiting pro-apoptotic Bad [66]. Akt 
activity is controlled by growth factors, which stimu- 
late membrane receptors and induce activation of Ras, 
transmitting signal to PI3K, which in turn phosphory- 
lates PIP2 into PIP3 [51,67] (PI3K can respond to growth 
factors also independently of Ras [68]). PIP3 enables 
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membrane localization of Akt, allowing for Akt activa- 
tion via phosphorylation at Thr308 and Ser473 by kinase 
PDK1 [69]. The anti-apoptotic Akt and its upstream 
regulators, such as GTPase Ras and kinase PI3K, are 
deregulated in a wide range of solid tumors and hema- 
tologic malignancies, hence the Akt pathway is consid- 
ered the key determinant of biological aggressiveness of 
these tumors and a major potential target for anticancer 
therapies [70,71]. 

Interestingly, phosphorylation of p53 at Ser46 enables 
it to activate expression of phosphatase PTEN [72,73], 
which prevents phosphorylation of Akt by dephosphory- 
lating PIP3 to PIP2. Only if activated, Akt mediates phos- 
phorylation of the p53 primary inhibitor, Mdm2, allowing 
it to localize to the nucleus and prime p53 for degradation 
[74,75]. These interactions intertwine tightly signaling of 
pro-apoptotic p53 and anti-apoptotic Akt. 

Apoptotic models 

Here we review mathematical models of the apoptotic 
pathway, which are relevant to our study. 

Stucki and Simon [76] have focused on inhibitors of 
apoptosis (IAPs) that are able to bind active caspases lead- 
ing to their degradation in the proteasome. They proposed 
a simple mathematical model, describing the molecular 
interactions between Smac, Smac deactivators, IAPs, and 
caspase-3, and derive the requirements for either induc- 
tion or prevention of apoptosis, which is initiated when 
the level of caspase-3 exceeds a given threshold. Further, 
Bagci et al. [77] described a mathematical representation 
of mitochondria-dependent apoptosis, in which kinetic 
cooperativity in the formation of the apoptosome (which 
consists of 7 cytochrome c-Apaf dimers) is a key ele- 
ment ensuring bistability in survival-or-death decisions. 
They examined the influence of Bax and Bcl-2 synthesis 
and degradation rates, as well as the number of mito- 
chondrial permeability transition pores (MPTPs) on the 
cell response to apoptotic stimuli. The model predicts 
that above some critical Bax degradation rate the sys- 
tem is monostable and cells survive. Bistability arises 
for Bax degradation smaller than critical, while the cell 
fate depends on the initial level of caspase-3. When 
the number of MPTPs is large, the bistability vanishes 
and the apoptosis is initiated regardless of the initial 
condition. 

Wee et al. [78] analyzed the mutual antagonism between 
pro-apoptotic signals of p53 and pro-survival signals 
of Akt in response to DNA damage. The coupling 
between p53 and Akt involves p53-regulated PTEN which 
dephosphorylates PIP3 required for Akt phosphorylation. 
The considered apoptotic module comprises four pro- 
teins from the Bcl-2 family: pro-apoptotic Bax and Bad 
(transcription of which is triggered by p53) and anti- 
apoptotic Bcl-2 and Bc1-xl (binding, respectively, Bax and 



Bad). In the model the pro-survival action of Akt relies 
on Bad phosphorylation which primes it for degradation. 
The authors demonstrated that repeated oscillations of 
the p53 level lead to the depletion of Bcl-2 and Bc1-xl, 
which is considered as a marker of apoptosis. Antagonism 
between pro-apoptotic signals of p53 and pro-survival sig- 
nals of Akt can also lead to bistability, which, if DNA 
repair is not accomplished in time, allows for termina- 
tion of p53 oscillations [79-81], significant increase of 
the p53 level and eventual apoptosis [82,83]. In corre- 
spondence to the mentioned studies [78,82,83], Li et al. 
[84] constructed an integrated mathematical model that 
includes three modules of the p53 network: p53 core reg- 
ulation, p53-induced cell cycle arrest and p53-dependent 
apoptosis initiation. Analysis of the model reveals that 
different aspects of the nuclear p53 dynamic profile 
are used to differentially regulate the pro-survival and 
pro-apoptotic modules. The initiation of caspase activity 
happens only when both the nuclear and the mitochon- 
drial p53 levels are above certain thresholds. However, 
the switch to the apoptotic state is mainly triggered by 
the accumulation of the mitochondrial p53, which pro- 
ceeds during oscillations of nuclear p53. Further, Tian 
et al. [85] proposed a two step mechanism coordinat- 
ing (early) mitochondrial and (late) nuclear p53 activities. 
In the model, apoptosis may quickly follow the severe 
DNA damage through the mitochondrial p53 pathway, 
or, in the case of less severe damage, apoptosis may fol- 
low the cycle arrest phase if DNA repair has not been 
accomplished. 

In summary, the above models revealed that apop- 
tosis is possible based on the system bistability, which 
makes the apoptotic decision irreversible. Bagci et al. 
[77] showed that bistability may arise due to cooperativ- 
ity in the apoptosome formation and two positive feed- 
back loops mediated by caspase-3. In the model of Tian 
et al. [85] bistability arises due to cooperation of two 
positive feedbacks: one is the double-negative feedback 
loop between p21 and caspase-3, the other positive feed- 
back loop is between cytochrome c and caspase-3. In 
correspondence to the model proposed here, the study 
of Wee et al. [78] indicated that the apoptotic deci- 
sion is controlled by opposing pro- and anti-apoptotic 
signals mediated respectively by p53 and Akt. Finally, 
Schlatter et al. [86] presented an apoptotic model based 
on the time-resolved Boolean network with multi-value 
node logic, which allowed them to emulate typical apop- 
totic features. The model is analyzed with regard to its 
internal connectivity and crosstalks, with special atten- 
tion to feedback loops and delayed processes such as gene 
regulation. 

In this study we aim to address the question, whether 
the interaction network of Bcl-2 family proteins allows for 
the integration of apoptotic signals, and if yes, whether 
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this signal integration is analogous to Boolean logic gates. 
Boolean networks are used to model various regulatory 
pathways, including apoptotic ones [86], however, the 
analysis of correspondence between biochemical reac- 
tion kinetics (which can be approximated by systems of 
ODEs) and the Boolean approach is missing. Here, fol- 
lowing Wee et al. [78] we focus on two key pathways 
leading to apoptosis: one mediated by p53 and the other 
by Akt As already discussed, p53 activates expression 
of PTEN, thus the elevated level of p53 results in the 
decreased level of phosphorylated Akt. However, since 
Akt activity is also regulated by growth factors (Figure 1), 
in our model the p53 and Akt pathways will be consid- 
ered independently. We will analyze how the apoptotic 
signals are collected and processed before the apoptosis- 
or-survival decision is attained. In accordance with pre- 
vious studies, we will assume that the apoptotic switch 
is based on bistability arising due to positive feedback 
mediated by caspase-3, and nonlinearity. We will demon- 
strate that the topology of the mitochondrial apoptotic 
module allows for the integration of signals in the man- 
ner analogous to the logic gates AND or OR, depend- 
ing on the levels of pro-apoptotic Bad and pro-survival 
Bcl-x L . 

Methods 

The system of ODEs following the proposed model 
was solved using the ode2 3tb solver implemented in 
Matlab (The MathWorks Inc., Natick, MA, USA). The 
source code is available as the Additional file 1. 

Results 

Model 

Our model of the apoptotic decision module involves 
the following components: p53, Akt, scaffolding pro- 
tein 14-3-3, caspases, and the Bcl-2 family proteins. The 



evidence that p53 serves as a transcription factor for bax 
is well-established [22]. Contrary to Bak, Bax can interact 
with any pro-survival restrainer [36]. The competition for 
Bc1-xl was reported between Bax and Bad, with the 
latter having higher affinity [42]. For these reasons, pro- 
apoptotic effectors Bax and Bak are represented col- 
lectively by a single entity termed Bax, and Bc1-xl is 
selected to represent the group of restrainers. Neutraliza- 
tion of Bc1-xl is the primary pro-apoptotic function of Bad 
[87,88], and it was demonstrated that a single mutation in 
Bad is sufficient to disrupt Bad:Bcl-2 but not Bad:Bcl-XL 
binding [88]; Bad binds to restrainers stronger than e.g. 
Bid [40]. Thus, BH3-only proteins are represented by a 
single Bad. 

Active Bax, when freed from its antagonists, releases 
cytochrome c from the mitochondrial intermembrane 
space [28]. Released cytochrome c enables the formation 
of the apoptosome, which in turn initiates the down- 
stream caspase program [10]. In the model, we omit 
these steps and assume that free Bax induces cleavage 
of pro-caspases into active caspases. The major caspase 
involved in this step is caspase-9, which can cleave itself 
into the active form within the apoptosome, and then 
activates the main executor caspase-3. Caspase-3 cleaves 
Bcl-2 leading to the further release of cytochrome c 
creating positive feedback loop, which introduces bista- 
bility and makes the apoptotic decision irreversible 
[14,15,83]. In the model, for the sake of simplicity, we 
consider a single caspase species and we simplify the 
positive feedback loop to the caspase auto-activation 
mechanism. 

In the model (Figure 2), levels of phosphorylated Akt 
(named here Akt p ) and p53 in its killer form (i.e. p53 phos- 
phorylated at Serl5, Ser20 and Ser46; denoted p53kiii e r) 
will serve as inputs. Surviving cells are characterized by 
high level of Akt p and lack (or very low level) of p53kin e r. 
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Figure 1 Crosstalk of pro-apoptotic and pro-survival pathways. The intracellular pro-apoptotic signal is mediated through p53kiiier- Elevated 
p53kiiier activates expression of PTEN, which leads to the inactivation of pro-survival Akt. Activity of Akt is also regulated by extracellular growth 
factors. For the apoptotic module analysis, both input signals, p53kiii e r and Akt Ul are considered independently. Arrow-headed lines indicate protein 
transformation, circle-headed lines denote activation or transcriptional regulation, hammer-headed lines denote repression. 
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Figure 2 Mitochondrial apoptotic module. The levels of p53kiiier and Akt u are considered as pro-apoptotic input signals. Total level of Akt is 
assumed constant and, since the phosphorylated Akt is anti-apoptotic, the unphosphorylated Akt is considered as the pro-apoptotic input 
(associated with e.g. growth factor withdrawal). Phosphorylated Akt phosphorylates Bad, which then remains bound to Scaffold 14-3-3. When 
dephosphorylated, Bad binds to Bcl-x L , leading to the dissociation of Bax from its complexes with Bcl-x L . Free Bax (via the omitted process of 
cytochrome c release) turns pro-caspases into active caspases (module output) which induce apoptosis. Arrow-headed lines indicate protein 
transformations, circle-headed lines - activating (including transcription) or inactivating influence. Light red color indicates pro-apoptotic, and light 
blue - anti-apoptotic proteins. Symbol 0 is used to denote sink or source. 



In these cells, most of Bad remains in the phosphory- 
lated form (denoted Bad p ) bound to Scaffoldi4_3_3, while 
Bax is inhibited (sequestered) by BcI-xl [44]. Unphos- 
phorylated Bad (Bad u ) may bind to Bc1-xl, which limits 
the amount of Bc1-xl protein available for inhibiting Bax. 
The steady state is controlled by the total levels of Bad 
(Badtot) and Bc1-xl (BcI-xl tot)> which will be considered 
as parameters in the model. The fraction of BcI-xl bound 
to Bad increases with Bad to t rendering cells of high Bad 
level more prone to apoptosis. (Protein levels, written in 
italics, are given in molecules per cell; for disambigua- 
tion, levels of protein complexes will be written in curly 
braces.) 

We consider two pro-apoptotic stimuli, which may lead 
to the accumulation of free Bax and consequent apoptosis: 

1) increase of the level of p53kin er , which triggers Bax 
transcription and Bax protein accumulation; 

2) dephosphorylation of Akt, i.e. decrease of Akt p and 
corresponding increase of the level of 
unphosphorylated Akt (Akt u ). 

Akt dephosphorylation is followed by the dephosphory- 
lation of Bad and its release from Scaffoldi4_3_3. Dephos- 
phorylated Bad captures BcI-xl which eventually releases 
Bax. 

Important notice: We assume that the total level of 
Akt {Akttot = Aktp + Akt u ) remains constant. Under 



this assumption, the appearance of unphosphosphory- 
lated Akt (Akt u ) will be considered as a pro-apoptotic 
stimulus. The strength of both pro-apoptotic stimuli will 
be thus measured by p53\^[\\ er and Akt u . 

The mathematical representation of the apoptotic 
module consists of 11 ordinary differential equations 
for levels of molecular species or their complexes: 
Bax m RNA> Bax, BcI-xl, Bax:Bcl-XL, Bcl-XL:Bad u , Bad u , 
Badp, Bad p :Scaffoldi4_3_3, Scaffoldi4_3_3, pro-caspases and 
caspases. 



^53killer 



M 2 +^53 ki ii er 2 

— d\Bax m KNA> 
^Bax(t) = s^BaXmRNA + u\{Bax\Bcl-xi} 

— b\Bax-Bcl-XL — d^Bax, 
^Bcl-Xiit) = U2{Bcl-Xi:Bad u } + u\{Bax\Bcl-xi} 

+ piAkt p {Bcl-Xi :Bad u } 

— b2Bcl-XL-Bad u — b\Bax-Bcl-xi 
+ d2{Bax:Bcl-Xi}, 

^{Bax:Bcl-xi}(t) = b\Bax-Bcl-xi — u\{Bax\Bcl-xi\ 

— d2{Bax:Bcl-Xi}, 

Bad u }(t) = b2Bcl-XL-Bad u — U2{Bcl-xi: Bad u } 
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— p\ Aktp {BcI-xl • Bad u }, 
^Bad u (t) — U2{Bcl-xi : Bad u ] — b^Bcl-xi -Bad u 

— piAktp -Bad u + q\Bad p 
+ qi {Badp : Scaffold 14-3-3 }, 

^Badp(t) = us{Badpi Scaffoldi4,-3-3} 

— b^Badp- Scaffoldings 
+ piAktp-Bad u 

+piAktp{Bcl-XL:Bad u } —qiBadp, 
^{Badpi Scoff old m-3-?>}(t) = bsBadp'Scaffoldmss 

— us {Badp : Scaffoldings } 

— q\{Badp : Scaffoldm-3-3 }, 
^Scaffoldmssit) = us{Badpi Scaffoldmss} 

— bsBadp'Scaffoldinss 

+ q\\Badp : Scaffoldm-3-3 }> 
^Procasp(t) = S3 — a\Bax-Procasp 

— a^Casp 1 -Procasp 

— d^Procasp, 
^Casp(t) — aiBax-Procasp 

+ a^Casp 1 -Procasp 

— d^Casp. 

Bax m RNA transcription, induced by pSSkuier, is modeled 
according to the sigmoidal kinetics (Hill equation with coop- 
erativity coefficient equal 2) [89]. The cubic term in caspase 
activation introduces bistability and assures irreversibility of 
the apoptotic decision [15]. Model parameters are provided in 
Table 1. 

The caspase switch 

The last two equations exhibit bistability. In Figure 3 
we show the bifurcation diagram of Casp with Bax con- 
sidered as a bifurcation parameter. Casp undergoes the 
saddle-node bifurcation at Bax = Bax^a — 5000 
(with Casp ~ 1000). Below the bifurcation point there exist 
three steady states: one unstable and two stable corresponding 
to low (< 1000) and high (> 9.9 x 10 4 ) caspase levels. Above 
the bifurcation point there exists the unique stable steady 
state characterized by high Casp. The structure of the bifurca- 
tion diagram ensures that the apoptotic switch is irreversible, 
i.e. once the system switches to the apoptotic state, it may not 
switch back even if the level of free effectors, Bax, drops to 
zero. The apoptotic switch occurs when Bax exceeds Bax^a for 
sufficiently long time; as we will show later, short excursions 
of Bax over Bax D [f do not induce apoptosis. For the steady 



state analysis, however, the states with Bax > Bax^a will 
be interpreted as apoptotic. 

OR and AND logic gates 

We demonstrate that apoptosis can be controlled in a 
manner similar to logic gates OR and AND with inputs 
defined as ^53kiii e r and Akt n . In this section, we con- 
sider digital steady-state inputs: the signal from p53kin er 
is assumed to be logic 0 for ^53^1^ = 0, and is assumed 
logic 1 for ^53kiiier = ^53j^fg r = 2 x 10 5 , which is the 
highest level considered in the analysis. The signal from 
the Akt branch is interpreted as logic 0 for Akt u = 0 
(i.e. when all Akt is in the phosphorylated form) and is 
considered logic 1 when Akt u = Akt to t = 2 x 10 5 . 
Intermediate input values of p53\^\\ eY e (0, i^53^g r ) and 
Akt u e (0, Akttot) will be analyzed in a forthcoming 
section. 

We will analyze the steady states of the system for two 
levels of Bad in order to demonstrate that the gate OR is 
achieved for Bad to t = 2 x 10 5 , while gate AND is achieved 
for Bad^t = 0-6 x 10 5 ; here, for both gates we assume 
the same BcI-xl tot = 1 x 10 5 . Surviving cell steady state is 
associated with conditions in which both apoptotic signals 
are equal zero. In these cells (Figures 4A and 5A) most of 
Bad is in the phosphorylated form bound to Scaffoldi4_3_3. 
The remaining (unphosphorylated) Bad is bound to 
BcI-xl. Rest of BcI-xl sequesters Bax or remains free. Free 
BcI-xl may be considered as the anti-apoptotic buffer, 
which may potentially capture the excess of Bax. This 
buffer of free BcI-xl is small (about 10 4 molecules) for gate 
OR (Figure 4A) and much larger (about 6 x 10 4 molecules) 
for gate AND (Figure 5A). Size of the BcI-xl buffer deter- 
mines cell susceptibility to apoptosis: small buffer renders 
cells more prone to apoptosis. 

Accordingly, to trigger apoptosis, gate OR (Figure 4) 
requires any of two signals. The signal coming from p53 
(Figure 4B) causes that the total Bax level increases such 
that Baxtot > Bcl-xi to t — Bad u . As a result, not all Bax 
may be sequestered by BcI-xl and free Bax appears. The 
signal coming from Akt (whole Akt desphosphorylation, 
Figure 4C) causes dephosphorylation of Bad, which is 
released from Scaffoldi4_3_3 and captures the whole pool 
of BcI-xl. As a result, all Bax (2 x 10 4 molecules) is 
released. 

Gate AND arises when the total level of Bad is smaller 
and correspondingly the initial (surviving cell steady state) 
level of (free) BcI-xl is larger (Figure 5, Table 2). In this 
case, both pro-apoptotic signals are required to trigger 
apoptosis. Although the p53 signal increases the level of 
Bax, still Bax^t < Bcl-xitot — Bad U) and thus all Bax 
remains sequestered by BcI-xl (Figure 5B). The signal 
coming from the Akt branch leads to dephosphoryla- 
tion of Bad, which dissociates from Scaffoldi4_3_3 and 
binds to BcI-xl (Figure 5C). However, again Bax to t < 
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Table 1 Model parameters 


Symbol 


Value 


Unit 


Description 


S] 


io- 2 


mlcl x s _1 


Basal Bax mRNA synthesis rate 


52 


3 x IO" 2 


mlcl x s _1 


p53kiiierregulated Bax m RNA synthesis rate 


53 


2 x 10 1 


mlcl x s _1 


Pro-caspases synthesis rate 


54 


2 x 1(T 1 


s- 1 


Bax protein synthesis 


M 


10 5 


mlcl 


Michaelis-Menten coefficient for p53kiii e rregulated Bax transcription 


di 


10- 3 


s~ 1 


Bax mRNA degradation rate 


d 2 


10- 4 


s~ 1 


Bax degradation rate 


d 3 


2 x i(r 4 


s~ 1 


Pro-caspases and caspases degradation rate 




3 x IO" 5 


mlcr 1 s- 1 


Bax-Bcl-XL binding rate 


b 2 


3 x i(r 3 


mlcr 1 s- 1 


Bcl-XL-Bad u binding rate 


b 3 


3 x i(r 3 


mlcl -1 s _1 


Badp-Scaffold 14-3-3 binding rate 


U] 


1Q -4 


s~ 1 


Bax:Bcl-XL heterodimer unbinding rate 


u 2 


10- 4 


s~ 1 


Bcl-XL:Bad u heterodimer unbinding rate 


u 3 


1Q -4 


s~ 1 


Bad p :Scaffoldi4_3-3 heterodimer unbinding rate 


Pi 


3 x i(r 10 


s~ 1 


Bad u phosphorylation rate in Bcl-XL:Bad u heterodimer (by Akt p ) 


<7i 


3 x 10~ 5 


s~ 1 


Badp dephosphorylation rate 


0] 


2 x 1CT 10 


mlcl -1 s _1 


Pro-caspases activation rate (by Bax) 


02 


10- 12 


mlcr 2 s- 1 


Pro-caspases autoactivation rate 


Akt X ox 


2 x 10 5 


mlcl 


Constant pool of Akt,/\/ct to t = Akt u + Akt p 


/-^killer 


2 x 10 5 


mlcl 


Maximum allowed level of p53kiiier 



Abbreviation 'mlcl' denotes the number of molecules per cell. 



Alternatively, gate OR may be transformed into gate 
AND (denoted AND* to avoid confusion) by a signifi- 
cant increase of the level of pro-survival Bc1-xl (Figure 6, 
Table 2). When Bcl-x Ltot = 2.4 x 10 5 (with Bad tot = 
2 x 10 5 ), dephosphorylated Bad is not able to displace Bax 
from its complex with Bc1-xl: high abundance of Bc1-xl 
allows it to restrain whole Bad u and whole Bax at the same 
time. 

In the next section we summarize the above numerical 
analysis introducing approximate formulas defining OR 
and AND gates. 

Macro-parameters characterizing steady-state solutions 

An important property of the model is that, as long as 
steady-state solutions are considered, a change of any of 
its 20 parameters can be compensated by a proper change 
of other (conjugate) parameter(s). Therefore, steady-state 
concentrations of proteins depend on clusters of param- 
eters rather than on individual parameter values. We 
identified 7 such parameter clusters, which define 'macro- 
parameters' (Table 3). When considered pairwise, macro- 
parameters are mostly mutually independent. 

Based on the macro-parameters one can obtain approx- 
imate algebraic formulas which must be satisfied for 
a given gate type. First, let us notice that since all 
(non-dimensional) dissociation constants l/m^, l/m^, 



BcI-xl tot — Bad u and thus all Bax remains captured by 
BcI-xl- Only the combination of signals coming from 
p53 and Akt, which leads to the simultaneous increase 
of the total Bax level and partial sequestration of BcI-xl 
by Bad u , results in the release of Bax (since Baxtot > 
BcI-xl tot — Bad u ) and eventually triggers apoptosis 
(Figure 5D). 
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10000 



1000 



100 



10 



1 1 1 1 




APOPTOSIS : 







2000 



4000 6000 
Bax 



8000 



10000 



Figure 3 Bifurcation diagram for Casp vs. Bax. Saddle-node 
bifurcation point is {Bax^,Caspw) ~ (5000, 1000). The unstable 
steady state is marked by dashed line. Solid lines show high and low 
stable steady states corresponding to apoptosis and survival, 
respectively. (Note the logarithmic scale on the vertical axis). 
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Apoptotic logic gate OR (Bad tot = 2 x 10 5 , Bcl-x Ltot = 1 x 10 5 ) 



No DNA damage: 

P 53 killer = 0 

High GF stimulation: 
AkL = 0 



B 



DNA damage: 

P 53 killer = P 53 S?kr 

High GF stimulation: 
AkL = 0 



No DNA damage: 

P 53 killer = 0 

GF withdrawal: 

AkL — Aktirt 



D 



DNA damage: 

P 53 killer = P 53 idfler 

GF withdrawal: 

A.kt u — A.ktf n f 



200 000 



130 000 



100 000 
90 000 



70 000 



0 





w 



























200 000 



140 000 
130 000 



100 000 



70 000 



200 000 



120 000 
100 000 



BaxBcl-x L Bad u Bad p Scaffold 14 _ 3 
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o 



200 000 



170 000 



100 000 1 



Bax Bcl-x L Bad u Bad p Scaffold 14 _ 3 

{0,1} => APOPTOSIS 




Bax Bcl-x L Bad u Bad p Scaffold 14 . 3 

{1,1} => APOPTOSIS 



BaxBcl-x L Bad u Bad p ScafFold 14 . 3 . 

{1,0} => APOPTOSIS 

Figure 4 Steady-state concentrations for gate OR. The presence of a single sufficiently strong input signal, p53kiiier or Akt Ul triggers apoptosis. 
Heights of bars correspond to the steady-state levels of proteins. Bars sharing adjacent sides (marked with colons and color gradient) denote 
complexes. To keep the convention, some bars do not start from 0. Accordingly, in (A) there are about 2 x 1 0 4 molecules of Bax:Bcl-XL complex, 
about 1 x 10 4 molecules of free BcI-xl, about 7 x 10 4 molecules of Bcl-XL:Bad u complex, about 13 x 10 4 Bad p :Scaffold 14-3-3 complex, and about 
7 x 10 4 molecules of free Scaffold 14-3-3. (A) The steady state levels of proteins and complexes corresponding to p53kiiier = 0 and Akt u = 0 (surviving 
cell, in the presence of growth factors). In order to induce apoptosis, gate OR requires a single sufficiently strong signal from (B) p53kiiier 
(p53kiiier = 2 x 1 0 5 ) or (C) Akt u {Akt u =2x1 0 5 ). (D) Simultaneous presence of both input signals (p53kiiier = 2 x 1 0 5 , Akt u = 2 x 1 0 5 ) results in 
higher level of Bax than in the case of stimulation by a single signal. 



and l/ms are larger than the characteristic numbers of 
molecules, all allowed complexes are formed. The condi- 
tion ms implies that BcI-xl binds preferentially 
to Bad u and therefore free Bax is released only when 

Baxtot + Bad u > BcI-xl. 

When the level of free Bax exceeds Bax D [f, the cell sets 
in the apoptotic state; in the opposite case the cell sur- 
vives. The amount of Bad u is (1 — me)Bad to v In this way 
we obtain simple inequalities defining logic gates. These 
inequalities, given in Tables 4 and 5, show the connec- 
tion between ODEs and Boolean logic gates, and pose 
restrains on parameter values. 

Level of Bad (BcI-xl) controls the transition between the 
logic gate AND and OR (AND* and OR) 

In Figure 7 we investigate system responses to ^53kin er and 
Akt u both in the range (0, 2 x 10 5 ), assuming levels of 
Badtot and Bcl-xi tot adequate to the gate type (Table 2). 
The Bax = Bax D a = 5000 isoline (white line) separates 
inputs {/?53kiiier> Akt u } leading either to apoptosis or to sur- 
vival. In gate OR (Figure 7A), apoptosis is activated for 



a relatively weak stimulation. It is triggered when steady 
state inputs satisfy j?53kiii e r > 0.66 x 10 5 or Akt u > 0.65 x 
10 5 , but even weaker signals, when in cooperation, may 
result in apoptosis. In contrast, in gate AND (Figure 7B) 
apoptosis requires cooperation of two stronger signals; 
for either ^53 ki iier < 0.84 x 10 5 or Akt u < 1.37 x 10 5 
apoptosis cannot be initiated. In gate AND*, a higher 
Akt u than in gate AND is required to trigger apoptosis 
(Figure 7C). 

In Figure 8A we determined Bax = 5000 isolines in 
the (p53k[\\ er , A/:£ u )-plane for different levels of Bad to t- For 
{^53kiiie r) Akt u } "above" each isoline value the cell under- 
goes apoptosis, whenever its Badtot is greater or equal 
Badtot value for that isoline. Gate AND arises for Badtot ^ 
0.9 x 10 5 , while gate OR arises for Bad tot > 1.1 x 10 5 . For 
Badtot € (0.9 x 10 5 , 1.1 x 10 5 ) the full dephosphorylation 
of Akt (Akt u = Akttot) leads to apoptosis (for an arbitrary 
j?53kiiier) but the increase of j?53kiii e r to its highest assumed 
value (/?53£Jfk r ) does not suffice for apoptosis without the 
additional signal from Akt. 

The transition from OR to AND* gate can be analyzed 
analogously by varying Bcl-xi tot (Figure 8B, Table 2). 
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Apoptotic logic gate AND (Bad tot = 0.6 x 10 5 , Bcl-x Ltot = 1 x 10 5 ) 
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Bax Bcl-x L Bad u Bad p Scaffold 14 _ 3 _3 

{0,1} => SURVIVAL 
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Bax Bcl-x L Bad u Bad p Scaffold 14 _3_ 3 

{1,1} => APOPTOSIS 



Bax Bcl-x L Bad u Bad p Scaffold 14 _3_ 3 Bax Bcl-x L Bad u Bad p Scaffold 14 _3_ 3 

{0,0} => SURVIVAL {1,0} => SURVIVAL 

Figure 5 Steady-state concentrations for gate AND. Transition from OR to AND results from the decrease of 8c?d t ot-The simultaneous presence 
of both input signals, p53kiii e r and Akt Ul is necessary to induce apoptosis. Heights of bars correspond to the steady-state levels of proteins. Bars 
sharing adjacent sides (marked with colons and color gradient) denote complexes. To keep the convention, some bars do not start from 0. (A) The 
steady state levels of proteins and complexes corresponding to p53kiii e r = 0 and Akt u = 0 (surviving cell, in the presence of growth factors). (B, C) a 
single maximal pro-apoptotic signal neither from p53kiii er nor from Akt u suffices to trigger apoptosis. (D) Simultaneous presence of both signals 
triggers apoptosis. 



Apoptosis and survival in response to transient stimulation 

In this section we determine the minimum duration of 
the stimulation phase (i.e. the phase in which one or 
two input signals {p53k[\\ er , Akt u } are present) needed to 
trigger apoptosis. We assume that before and after the 
stimulation phase p53\^\\ er = 0 and Akt u = 0. The 
transient stimulation by pSS^iiier and/or Akt u may cause 
that the caspase level passes a threshold above which the 
apoptotic decision is irreversible. This threshold caspase 
level 2 x 10 3 ) is determined by the value of the unsta- 
ble steady state for Bax = 0, see the bifurcation diagram 
(Figure 3). 



Table 2 Apoptotic gate types resulting from considered 
Bcl-2 family protein levels 





Protein levels 




Gate type 


Baxtot 


BcI-xl tot 




dynamic 


1 x 10 5 


2 x 10 5 


OR 


dynamic 


1 x 10 5 


0.6 x 10 5 


AND 


dynamic 


2.4 x 1 0 5 


2 x 10 5 


AND* 



Levels of total BcI-xl and Bad are assumed constant while the level of total Bax is 
dynamically regulated. 



For the gate AND we consider stimulation during which 

Akt u = Akttot and /?53killer = ^ 53 kiUer' We l <now fr° m 
the analysis presented in Figure 8 A (point {1,1}) that 
such stimulation lasting sufficiently long leads to apopto- 
sis. Here we estimated that the minimal duration of the 
stimulation phase is 10.5 hours. In Figure 9A we show 
that 10 hour long stimulation is insufficient for trigger- 
ing apoptosis; despite Bax reaches 1.3 x 10 4 > Baxuf 
and Casp reaches 1.6 x 10 3 > Caspbtf, the system returns 
to the initial steady state of the low caspase level. In 
Figure 9B we show that after 11 hour long stimulation 
the caspase level passes the unstable steady state lead- 
ing to the irreversible transition after which the caspases 
level settles at the high stable steady state. It is worth 
noting that after the stimulation Bax returns to its initial 
steady state. In the case of gate AND*, the minimal dura- 
tion of the stimulation phase (for the same stimulation) is 
20.5 hours. 

For gate OR we consider three particular modes of stim- 
ulation (Figure 10). In the first case we estimated that the 
minimum duration of the stimulation phase for the sole 
/*53kiiier signal ({^53 ki ii e r = /*53gg£ r , Akt u = 0}, point {1,0} 
in Figure 8A) is 3.0 hours. In the case of Akt-only signaling 
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Apoptotic logic gate AND* (Bad tot = 2 x 10 5 , Bcl-x Ltot = 2.4 x 10 5 ) 
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P 53 killer = 0 

High GF stimulation: 
Akt u = 0 
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{0,0} => SURVIVAL {1,0} => SURVIVAL {0,1} => SURVIVAL {1,1} => APOPTOSIS 

Figure 6 Steady-state concentrations for gate AND*. Transition from OR to AND* results from the increase of £c/-^ to t-The simultaneous 
presence of both input signals, p53kiii er and Akt Ul is necessary to induce apoptosis. Heights of bars correspond to the steady state levels of proteins. 
Bars sharing adjacent sides (marked with colons and color gradient) denote complexes. To keep the convention, some bars do not start from 0. (A) 
The steady state levels of proteins and complexes corresponding to p53knier = 0 and Akt u = 0 (surviving cell, in the presence of growth factors). (B, 
C) A single maximal pro-apoptotic signal neither from p53kiii er nor from Akt u suffices to trigger apoptosis. (D) Simultaneous presence of both 
signals triggers apoptosis. 



({/^kiiier = 0, Akt u = Akt tot }, point {0,1} in Figure 8 A) 
the critical duration of the stimulation phase is 2.8 hours. 
A slightly shorter stimulation (approximately 2.7 hours) 
is required when both signals are present simultaneously, 
even when their amplitudes are twice smaller than in the 



previous case ({p53\^[\\ er 

US 



max a i f 
killer' A/Cr u 



\Akt u 



point in Figure 8A). In Figure 10 we show trajecto- 



ries corresponding to the transition to apoptosis following 
3 hours long stimulation in each of three cases discussed 
above. 



Table 3 Macro-parameters governing steady-state solution levels 

Macro-parameter definition Value Unit 



Description 



my ■ 
m 2 :■■ 



51 5 4 
djd2 



" V 2m2+(p53 »V 



2 x 10 4 
6.8 x 1 0 4 



mid 
mlcl 



m 3 := 



0.3 



rnlcr 



Bqx-BcI-xl affinity 



/7?4 := 



m 5 := 



— k 



30 
30 



mlcl- 
mlcr 



Bcl-XL-Bad u affinity 

Bo dp-Scoffo Id] 4-3-3 affinity 



m 6 := 



m 7 := 



P] Akt tot 
q]+P] Akhot 
d 3 S-a 2 S 2 (^-S)~ 



where S = Cosp + Procosp 



= 5000 



mlcl 



max. fraction of phosphorylated E 
Box b]f (at Cosp = 1 000) 



Abbreviation 'mlcl' denotes the number of molecules per cell. 
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Table 4 Algebraic conditions for outputs of the OR gate 



Gate OR 


Inputs 


Output 


Condition for the output 


{0,0} 


Survival 


nri] + (1 — me)Bodt 0 t < Bcl-xi + m 7 


{1,0} 


Apoptosis 


nri2 + 0 — m^Badtot > BcI-xl + m 7 


{0,1} 


Apoptosis 


rri] + Bad ^ > Bcl-x L + m 7 


{1,1} 


Apoptosis 


m 2 + Badtot > Bcl-x L + m 7 



Input in the format {/J} denotes p53kiii er = / x2 x '\0 5 ,Akt u =j x 2 x 10 5 . It is 
assumed that m 3 m 4 . 



Discussion 

Robust response to genotoxic stress is crucial for cell 
development and enables cells to preserve their genome 
integrity and suppress oncogenic transformation [2]. We 
proposed the ODE model of the apoptotic module, 
allowing for unambiguous cell death or survival deci- 
sions. The model was constructed on the basis of our 
current understanding of interactions of Bcl-2 family 
members in regulation of cell fate in response to stim- 
uli mediated by two important upstream pathways of p53 
and Akt In the wider context, apoptosis is initiated in 
response to various stress stimuli, which often crosstalk 
with anti-apoptotic signals [90]; it is thus important to 
recognize, also from the mathematical point of view, 
how these signals might be integrated into survive-or-die 
decisions. 

The Boolean modeling, which typically utilizes the qual- 
itative knowledge about regulatory pathways, is often 
considered as a preliminary step to continuous models 
which require more detailed (kinetic) data, but allows also 
for a more precise description. Wittman and colleagues 
[91] introduced a systematic way of transforming Boolean 
models into ODE models by employing multivariate poly- 
nomial interpolation. Here, we show the reverse cor- 
respondence, demonstrating that the ODE model may 
behave akin to logic gates, in which YES or NO responses 
are associated with the change of relative abundances 
of proteins from pro- and anti-apoptotic subgroups of 
the Bcl-2 family, allowing for the release of the more 



Table 5 Algebraic conditions for outputs of the AND and 
AND* gates 

Gates AND and AND* 



Inputs 


Output 


Condition for the output 


{0,0} 


Survival 


m] + (1 — me)Badtot < BcI-xl + m 7 


{1,0} 


Survival 


m 2 + (1 - m 6 )Bad lo t < Bcl-x L + m 7 


{0,1} 


Survival 


my + Badtot < Bcl-x L + m 7 


{1,1} 


Apoptosis 


m 2 + Badtot > Bcl-x L + m 7 



Input in the format {i,j} denotes p53 killer = / x 2 x '\0 5 ,Akt u =j x 2 x 10 5 . It is 
assumed that rri3 1714. 



abundant protein. Interestingly, we found that our logic 
gates are reconfigurable: levels of pro-survival and pro- 
apoptotic proteins (Table 2) together with seven macro- 
parameters (assumed constant in the model; Table 3) 
determine gate type, AND or OR (by relations given 
in Tables 4 and 5). Another example of ^configura- 
bility was recently demonstrated by Goni-Moreno and 
Amos [92]. 

The considered mitochondrial apoptotic module inte- 
grates pro-apoptotic signals from p53kin er and pro- 
survival signals from Akt. The final output of the 
module is the level of caspase-3, which (in the model) 
is activated by Bax (or Bak), but also auto-catalytically, 
which renders the apoptotic decision irreversible. After 
the caspase level surpassed unstable steady state, the 
presence of Bax is no longer needed, since caspase 
autoactivation suffices for their further build-up and 
execution of apoptosis. Caspase activation requires 
that free Bax surpasses a threshold of about 5000 
molecules per cell, and remains over this level for 
sufficiently long time. The regulation of Bcl-2 fam- 
ily module can be summarized as follows: in resting 
cells, Bax (chosen as a representative for pro-apoptotic 
multidomain effectors) remains in the inactive form 
bound to BcI-xl (chosen as a representative for anti- 
apoptotic restrainers), while most of Bad (chosen as 
a representative for BH3-only pro-apoptotic proteins) 
is in the phosphorylated form bound to Scaffoldi4_3_3. 
The signal coming from p53 leads to the accumula- 
tion of Bax, while the signal coming from Akt (Akt 
dephosphorylation) results in Bad dephosphorylation. 
Dephosphorylated Bad may release Bax from BcI-xl. In 
this way, both mechanisms contribute to the appearance 
of free Bax. 

Cells characterized by high Bad level or/and low level 
of BcI-xl, have relatively small reservoir of free BcI-xl 
(required to inhibit/sequester Bax) and thus are more 
prone to apoptosis. In these cells apoptosis follows from 
p53kiiier accumulation or Akt dephosphorylation (gate 
OR); a BH3-domain mimetic, ABT-737, which mirrors 
binding capacities of Bad and engages pro-survival pro- 
teins (mostly Bcl-2), was shown to induce Bax/Bak- 
dependent killing [93]. Cells characterized by low Bad 
level or/and high level of BcI-xl have much larger reser- 
voir of free BcI-xl and in these cells apoptosis requires 
both signals simultaneously (gate AND); this is consistent 
with experimental results showing that the overexpres- 
sion of BcI-xl blocks apoptosis [94]. We demonstrated 
that transition between AND and OR gates results from 
either increase of Bad level or decrease of BcI-xl level, 
which confirms that levels of these proteins are important 
in regulating cell sensitivity to apoptosis. 

It is known that Bad modifications, which interfere with 
Bad phosphorylation, can make cells more or less sensitive 
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Figure 7 Bax in the input signals (p53kin er , ^4At u )-plane for three considered logic gates. (A) gate OR, (B) gate AND, (C) gate AND*. For all 
gates the black curve is the apoptotic isoline Box = Box^f = 5000. For values of Akt u and p53kiiier above the Box = Box^f isoline cell undergoes 
apoptosis, for values of Akt u and p53kiii er below this curve cell survives. 



to apoptosis. Bad phosphorylation is limited by PRMT1- 
mediated methylation of two Bad arginine residues (Arg94 
and Arg96), which prevents Akt-mediated phosphoryla- 
tion of Bad at Ser99. Respectively, decreased methyla- 
tion of Bad increases the fraction of phosphorylated Bad 



leading to its enhanced sequestration to Scaffoldi4_3_3, 
decreased caspase activity, and consequently enhanced 
cell viability [95,96]. 

In the therapeutic context, our findings suggest that in 
some (cancerous) cell lines or cell mutants (characterized 



Gate OR « Badtot > Gate AND 
Bax = 5000, Bcl-x L tot = 1 x 10 5 
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Figure 8 Transitions between gate types. (A) transition from gate OR to gate AND, (B) transition from gate OR to AND*. Apoptotic {Box = Box^f = 
5000) isolines are plotted in the input signals (p53y\\\ er , /\/ct u )-plane. Black dots in (A) denote {p53kiii e r, Akt u } pairs of input signals analyzed in Figures 9 
and 10. 
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Gate AND (Bad tot = 0.6 xlO 5 , Bcl-x Ltot = 1 x 10 5 ) 
{p53 kil i er = p53Sfier> Akt u = Akt tot } = {l,l} 5 active phase duration: 10 h 
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Gate OR (Bad tot = 2 xlO 5 , Bcl-x Ltot = 1 x 10 5 ), 
{p53 killer = p53gj£ r , Akt u = 0} = {1,0} 
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by the inhibition of pro-survival Akt pathway in order 
to be effective. In turn, it suggests that cells character- 
ized by high Bad level, or low Bc1-xl level are very prone 
to apoptosis, which can follow even from growth factor 
withdrawal. 

One could expect that levels or activities of proteins 
represented by Bad and Bc1-xl, which are assumed con- 
stant in the model, are also regulated in response to 
pro-apoptotic or pro-survival cues. Inclusion of these 
effects would expand the regulatory network and allow 
for the integration of a larger number of signals [86]. In 
order to improve the resolution of the presented apoptotic 
model, one should take into account individual charac- 
ters of proteins from the group of effectors, restrainers 
and, most importantly, BH3-only proteins, and in this way 
cover additional effects such as the induction of expres- 
sion of some BH3-only proteins by p53 [97,98], diverse 
specificities of BH3-only proteins with respect to vari- 
ous restrainers, and plausible requirement of the direct 
activation of effectors by some BH3-only proteins, which 
allows effectors to nucleate, oligomerize, and then release 
cytochrome c. 

Conclusion 

We demonstrated that the mitochondrial apoptotic mod- 
ule may process signals in the way similar to logic gates OR 
and AND, which suggests a mechanism of the integration 
of apoptotic and survival signals before the death- 
or-survival decision is reached. The correspondence 
between ODEs and Boolean logic arises due to the high 
affinity, competitive binding of two subfamilies of pro- 
apoptotic proteins, multidomain effectors (represented 
by Bax) and BH3-only proteins (represented by Bad) to 
anti-apoptotic restrainers (represented by Bc1-xl). Such 
regulation is based on stoichiometry of interacting pro- 
teins, and introduces dynamical thresholds for all protein 
levels, which, when surpassed, change output from NO to 
YES. Specifically, apoptosis arises when the level of anti- 
apoptotic restrainers is surpassed by the combined levels 
of pro-apoptotic unphosphorylated BH3-only proteins 
(controlled by Akt) and pro-apoptotic effectors 
(controlled by p53). We demonstrated that the transi- 
tion between OR and AND gates is accomplished by the 
change of a level of single component (either decrease 
of pro-apoptotic protein Bad or increase of pro-survival 
Bc1-xl), without any modification of topology of the 
network or even kinetic rate parameters. 



by a low Bad level or a high Bc1-xl level) apoptosis 
can result only from the simultaneous presence of both 
pro-apoptotic signals, i.e. elevated p53kin er level and Akt 
dephosphorylation (growth factor withdrawal). Radio- 
therapy against these cells which leads to the increase of 
p53kin er level (via DNA damage) should be accompanied 
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